Self-interaction spin effects in inspiralling compact binaries 



in 
o 
o 

(N 



> 
00 

m 
in 
^1- 
o 
in 
o 



Of 
i 1 

o 

CO 

c3 



X 



Balazs Mikoczi 1 , Matyas Vasiith 2 , Laszlo A. Gergely 1 
1 Departments of Theoretical and Experimental Physics, University of Szeged, Szeged 6720, Hungary and 
2 KFKI Research Institute for Particle and Nuclear Physics, Budapest 114, P.O.Box 49, H-1525 Hungary 

(Dated: February 2, 2008) 

Gravitational radiation drives compact binaries through an inspiral phase towards a final coales- 
cence. For binaries with spin, mass quadrupole and magnetic dipole moments, various contributions 
add to this process, which is characterized by the rate of increase df /dt of the gravitational wave 
frequency and the accumulated number N of gravitational wave cycles. We present here all contri- 
butions to df /dt and M up to the second post-Newtonian order. Among them we give for the first 
time the contributions due to the self-interaction of individual spins. These are shown to be com- 
mensurable with the proper spin-spin contributions for the recently discovered J0737-3039 double 
pulsar, and argued to represent the first corrections to the radiation reaction in the Lense-Thirring 
approach. 



I. INTRODUCTION 

Neutron-stars and black-holes forming compact binary 
systems emit gravitational radiation. The frequency 
range of this radiation is expected to be in the sensitivity 
range of the Earth-based interferometric detectors LIGO 
0, VIRGO @, GEO 0, TAMA H and currently ob- 
servation is under way to capture such signals. Recently, 
a method of setting upper limits on inspiral event rates 
for binary neutron stars using interferometer data were 
established from the first scientific run of LIGO and 
limits on gravitational wave emission from selected pul- 
sars using LIGO data were published 6]. Also analysis 
methods were specified for inspiral signals from binaries 
with 3-20 solar masses Gravitational radiation from 
these compact binaries also is expected to be detected by 
the Laser Interferometer Space Antenna (LISA) ||, [Ij. 
Parameters of the spinning compact binaries can be es- 
timated and alternative theories of gravity can be tested 
[Io| from these measurements. 

The final cataclysmic coalescence of such compact bi- 
naries is preceded by a milder inspiral phase, for which 
the post-Newtonian (PN) approach provides a reliable 
description. This description is generally considered valid 
until the system reaches the innermost stable circular or- 
bit (ISCO). For neutron star binaries with realistic equa- 
tion of state the gravitational wave frequency at ISCO is 
between 800 - 1230 Hz, according to [ill E3, O The 
upper limit of senstivity for the LIGO detector is ~ 1000 
Hz. The simple estimate shows that for a neutron star of 
1.4 solar masses the ISCO (at three Schwarzschild radii) 
is of about 2.5 times the characteristic radius of 10 km of 
the neutron star. A more sophisticated argument, tak- 
ing into account various equations of state for the neu- 
tron stars [3] shows risco/ r NS € (1.1,2.2). Buonanno, 
Chen and Vallisneri have discussed the problem of the 
failure of the PN expansion during the last stages of in- 
spiral (the intermediate binary black hole, IBBH problem 
|15|) both for non-spinning and for spinning black hole 
binaries on quasicircular orbits [0| UjJ • Reliable results 
can be achieved by a frequency cutoff, e.g. stopping the 
integration at the minimum of the energy as function of 



orbital frequency, at ICO [3 (called MECO in H3). In 
even latter stages of the inspiral, tidal torques become 
important |l9L |20j. Then each component can be rep- 
resented by a Schwarzschild metric corrected with the 
influence of the Weyl curvature generated by the other 
component 0, [2^ • The corrected metric was computed 
to second order in the inverse of the Weyl curvature ra- 
dius in [23L |24| and to third order by Poisson [2^ . 

In the post-Newtonian regime, the equations of motion 
were given to 3.5 PN order accuracy, with the inclusion 
of spin-orbit (SO) effects and their first PN correction 
in ptl . Spin-spin (SS) [27], quadrupole- monopole (QM) 
[2||and magnetic dipole - magnetic dipole (DD) contri- 
bution [2^| to the accelerations were also discussed. The 
backreaction on the orbit of the gravitational waves es- 
caping compact binaries can be characterized by the rate 
of radiative change of otherwise conserved quantities, like 
the secular energy loss (dE/dt) and secular angular mo- 
mentum loss (dJ/dt). These were computed over the 
years with great accuracy. The leading order contri- 
bution was given by J3fj|. |3lj . its first post-Newtonian 
(PN) correction in |32|. while the second post-Newtonian 
(2PN) correction in [33] (a correction is expected, see 
[p). 

Taking into account the interaction of the spins Si with 
the orbit, the total angular momentum J = L + Si + S2 
is still conserved, however the orbital angular momentum 
L is not, due to spin precessions. Still, its magnitude L 
remains a constant of motion, as consequence of the spe- 
cific functional form of the precession equation |35| . The 
radiative changes (dE/dt) and (dL/dt) characterize the 
backreaction on the radial part of the motion. They were 
computed in |35| and [36j (these approaches rely on dif- 
ferent spin supplementary conditions and distinct aver- 
aging techniques over a radial orbit - the results, however 
agree). The secular losses (dE/dt) and (dL/dt) depend 
on the relative angles k.; of the spins and orbital angular 
momentum. In order to have a closed system of differ- 
ential equations their radiative evolution (together with 
the radiative evolution of the ang le 7, subtended by the 
two spins) has been given in [35j |. For compact binaries 
the SO contribution occurs at e 3 / 2 order, e representing 
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one post-Newtonian order. 

The SS contribution to (dE/dt) and (dL/dt) and the 
subsequent evolution of the angles Ki and 7 were derived 
in j33, [3^|. Here a description in terms of the magni- 
tude of an (angular averaged) orbital angular momentum 
L was introduced, this convenient method being kept in 
the discussion of the QM contribution (3^ and the DD 
contribution |40j. (The magnitude L of the orbital an- 
gular momentum is not conserved due to the spin pre- 
cessions caused by spin-spin, quadrupole-monopole and 
magnetic dipole - magnetic dipole interactions.) The SS, 
QM and DD contributions all arise at e 2 order. 

Remarkably, the SS contributions in (dE/dt) and 
(dL/dt) given in [37j contained not only interaction 
terms between the two spins, but self-interaction spin 
terms (SS-self) as well. These arose from the terms 

proportional to JgQ 1 (&n) J so* 1 ( a Jv) in dE/dt and 
e ijkj(m ( aN )jW kl ( aw ) i n dJ 1 /dt (here J<0 1 denotes 
the n derivative of the spin-orbit contribution of the ve- 
locity quadrupole moment evaluated with the Newtonian 
acceleration a^v; there is summation understood over the 
repeated indices j,k,l). These terms are of order e 2 as 
well [H|. 

Although the SS-self contributions to the closed sys- 
tem of differential equations (d (E, L, K i: 7) /dt) has been 
given in [^3], H^, these terms are notoriously missing 
from the bookkeeping of different contributions to vari- 
ous results concerning inspiralling compact binaries. In 
this paper we compute the SS-self contributions to the 
rate of increase of the gravitational wave frequency / and 
to the accumulated number of gravitational wave cycles 
TV. For completeness we also enlist all other contribu- 
tions to df/dt and M to order e 2 (the PN, SO, SS, QM, 
DD, 2PN and tail contributions). 

On the long run due to the emission of gravitational 
waves the orbit tends to circularize [3(| ■ Therefore we 
consider quasicircular orbits, for which the gravitational 
wave frequency is twice the orbital freauencv|42j|. In 
Sec. II first we evaluate the rate of increase of /. This 
is given by the rate of change of the orbital angular fre- 
quency bj — 7r/ under radiation reaction: 



were given before in |3J, |33j, |35|-|40| and |43j. We have 
grouped together the contributions occuring at the same 
post-Newtonian order. 

Eq. Q being an ordinary differential equation in u 
can be immediately integrated. Then we compute the 
accumulated number of gravitational wave cycles TV by 
a second integration. Finally we evaluate these contri- 
butions to Af for some famous or typical compact binary 
systems discussed in the literature, and we compare the 
magnitude of the SS-self contribution with other contri- 
butions. This is done in the Concluding Remarks. 

In the following we define some of the notations em- 
ployed in the paper. 

The magnitude and direction of Si are denoted as 
Si and Si. Then 7 = cos _1 (Si • S2). The magni- 
tude and direction of the magnetic dipole moments are 
denoted as di and di. They subtend the angle A = 
cos _1 (di • CI2) with each other. In a coordinate systems 
K, with the axes (c, L x c, L), where c is the unit vector 
in the J x L direction, the polar angles Ki and ipi of the 
spins are defined as Si =(sin Ki cos ibj , sin Ki sin ipi, cos Ki). 
We also have introduced in [4(| coordinate systems 
JC 1 with the axes (bj, Si x bj, Si), where bj are the 
unit vectors in the S; x L directions, respectively. 
The polar angles on and /3j of the the magnetic 
dipole moments di in the coordinate systems JC 1 are 
di = (sin o>i cos j3i , sin on sin (3% , cos on ) . Quadrupolar con- 
tributions are characterized by the parameters Pi = 
Qi/mim 2 , where Qi is the quadrupole-moment scalar (2^ 
of the i th axially symmetric binary component with sym- 
metry axis Si. (In the Newtonian limit Qi = O,^ — Oj = 
— Si (Oi/O^ — 1) /Vli for each axisymmetric body charac- 
terized by the principal moments of inertia (O^, Oj, OJ) 
and angular velocity Sli — Si/O^.) 



Throughout the paper we use G 



1. 



dio 
~dt 



dE 
diu 



dE 
~dt 



(1) 



where the expression dE/du can be found by differen- 
tiating E = E(lo). We have verified that the circular 
orbit limit of the instantaneous energy loss (dE / dt) clrc 
is the same as the circular orbit limit of the secular en- 
ergy loss (dE / dt) clrc , the latter however is much simpler 
to compute. Here the various contributions to the secular 
energy loss 
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II. FREQUENCY EVOLUTION AND THE 
ACCUMULATED NUMBER OF CYCLES 



For circular orbits f = r = and v 



3 _ 



holds. Con- 



sequently (moj) 2 ^ is of order e. The radial projection 
of the acceleration defines the orbital angular velocity as 
r • a = — r 2 uj 2 [27j • From the explicit form of the acceler- 
ation (with various contributions given in |27| . [2§| and 
|29j) we then find lu = u> (r) and subsequently r = r (uj) 
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where 



"(w) = m(m^)- 2 / 3 Jl-^(mcc0 2/3 



f = ""SiSa + VSS-self + &QM + &DD 



(8) 



■ X E H+ 3! ? ^COSK, 

s \ 777/ / mf 



\4/3 



19 n 



+ - j (COS 7 — dCOSKi COS K2) 



2r/m 

2 



- (3 cos 2 - l) 



did 2 A 
2?7m 4 



(3) 



In the magnetic dipolar contribution we have introduced 
the shorthand notation: 



Aq = 2 cos A + 3(pi<72 - P2&1) sin(^2 - 

-3(pip 2 + CTiCr 2 ) COs(t/) 2 - ^l) , 



(4) 



The quantities 0, a SlS2 , (?ss-self, °~QM and a DD are the 
spin-orbit, spin-spin, self-interaction spin, quadrupole- 
monopole and magnetic dipole-dipole parameters, re- 
spectively: 



1 Si ( m 2 



' 9 =T^E^ H3^ + 75r, cos K4 , (9a) 



a SiS 2 



SiS: 



4877m 



-(-247cos7 + 721coskiCosk 2 ) , (9b) 



OSS-self = 7^-^ V ( — ) (6 + Sin 2 Ki) , (9c) 
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oqm =-7^]?! (3 cos 2 Ki - l) , (9d) 



with 



pi = sin aj cos /3j , 

(7, = cos a, sin + sin on sin /3, cos K4 



(5) 



Combining Eq. © with the expression of the velocity on 
circular orbits v = ru> in the energy integral E — E {r, v) 
we obtain E = E (uj): 

m = -i/i(mc) 2 / 3 (l-i(3+|)(m^) 2 / 3 
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The radiative orbital angular frequency evolution for the 
above defined circular orbits is deduced from Eqs. 0J, 
(J2J) and 10, with the various contributions given as 



o~dd — jdid2Ao , 



(9e) 



The N, PN, SO, SS, 2PN and tail contributions to Eq. 
Q were verified to agree with those given in 1441 . the 
QM with those in [2^] and the DD with those in |29j| , re- 
spectively. We emphasize the SS-self contribution, given 
for the first time. 

Eq. (JJJ is an ordinary differential equation in uj, pro- 
vided all angular variables can be considered constants. 
To linear order in the perturbations this is true. This 
is because all angles appear only in the perturbative cor- 
rections, therefore they are given with sufficient accuracy 
to e order, in which order they are constants [3{|. Thus 
the time-evolution of the orbital angular frequency u> for 
circular orbits can be easily deduced by an integration 
over time: 
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where the dimensionless time variable t = r)(t c — t)/5m is 
related to the time (t c — t) left until the final coalescence. 
A further integration gives the accumulated orbital phase 
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u)(r)dT 



(11) 
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From here the accumulated number of gravitational wave 
cycles emerges as 



Af 



= J_ J r 5/s , /3715 55 \ 3/8 
7T 7T77 1 V 8064 96 7 

+ 3 r i/4 | ( 9275495 284875 

4^4 / ^14450688 258048 ' 



2048 ' 64 i 
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The tail and 2PN contributions agree with those given in 
HI], the rest of the terms with Eq. (4.16) of l27l. W e also 
note that the numerical values in Table II of 10] could 
be recovered from our Eq. 112(1. 

In Tabledwe enlist all contributions to Af in terms of (3 
and a, evaluated numerically from our Eq. 1)12(1 for sev- 
eral compact binaries. These are two well-known binary 
neutron star systems, the newly discovered double pulsar 



J0737-3039 (fig, fig, and for a recent review see fig), 
the Hulse- Taylor pulsar B1913+16, one neutron star - 
stellar mass black hole binary, one binary composed of 
two stellar mass black holes |49|. finally two examples of 
galactic black hole - galactic black hole binaries [lCj . 



III. CONCLUDING REMARKS 



We have presented the complete set of contributions 
up to e 2 order (PN, SO, SS, QM, DD, tail and 2PN) to 
the evolution of gravitational wave frequency and to the 
accumulated number of gravitational wave cycles, with 
the inclusion of the previously unknown self-interaction 
spin terms. These results add to the already derived [35| , 
[13]- closed system of first order differential equations 
governing the secular evolution of radiating compact bi- 
naries and represent an important step towards a com- 
plete characterization of the orbital evolution. 



TABLE I: Number of accumulated gravitational wave cycles, 
either for the frequency domain of the detectors or up to the 
innermost stable circular orbit. The first four sets of data 
refer to the LIGO/ VIRGO type detectors, while the last two 
columns to the LISA bandwidth. 





J0737-3039 


B1913+16 


BDIWW1 


BDIWW2 


BBW1 


BBW2 


PN Order 


1.337M Q 


1.387M 


1.4M@ 


1OM 


1O 4 M 


1O 7 M 




1.250M q 


1.441M Q 


1OM 


10M Q 


1O 5 M 


1O 7 M 


MHz) 


10 


10 


10 


10 


4.199 x 10~ 4 


1.073 x 10~ 5 


ffin(Hz) 


1000 


1000 


360 


190 


3.997 x 10~ 2 


2.199 x 10~ 4 


N 


18310 


15772.1 


3580 


600 


21058 


535 


PN 


475.8 


435 


212 


59 


677 


55 


SO 


17.5/3 


16.5/3 


14/? 


4/3 


36/3 


4/3 


SS - self, SS, QM, DD 


-2.1(7 


-2.1a 


-3a 


— a 


-5a 


— a 


Tail 


-208 


-206 


-180 


-51 


-450 


-48 


2PN 


9.8 


9.5 


10 


4 


18 


4 



In what follows we would like to comment on the im- 
portance of the self-interaction spin contributions. For 
this, we evaluate the spin parameters <JS\S? and crss-self 
for the case of the recently discovered double pulsar 
J0737-3039. The neutron stars in this double pulsar have 
average radii 15 km, masses of 1.337 Mq and 1.25 M©, 
pulse periods of 22.7 ms and 2773.5 ms. According to the 
Jenet-Ransom model [5p| the angle k\ (in their notation 
A) has two possible values: Ki = 167° ± 10° (Solution 1), 
and Ki = 90° ± 10° (Solution 2). The other angle k 2 in 
principle can be determined by solving numerically the 



spin precession equations. This was done for black hole 
binaries in (SJ. However according to jEJ and |5^] it is 
likely that wind-torques from the energetically dominant 
component have driven the spin axis of the other com- 
ponent to align with the direction of L, causing «2 = 0. 
Therefore 7 = m. We give the estimates for the spin 
parameters in Table ITU 

Two remarkable conclusions stem out from Table [H] 
First, let us consider Solution 1. Though four orders 
of magnitude smaller than the spin-orbit effects, already 
considered in |54j . self- interaction spin contributions are 
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TABLE II: The spin parameters for the two solutions of the 
Jenet-Ransom model representing the binary pulsar J0737- 
3039. 



Spin parameter (order) 


Solution 1 


Solution 2 


P 


-0.166 


0.001 


°s lS2 (10- 4 ) 


-0.372 





O'SS-self (10~ 4 ) 


0.298 


0.345 



found to be comparable with the proper spin-spin con- 
tributions. This happens in spite of the fact that the 
numerical factors in ass-self are two orders of magni- 
tude smaller than those in <js 1 s 2 j ci - Eqs. (|9b|) and (|9"cj) . 
The commensurability of the proper and self SS contri- 
butions is due to the fact that one of the spins is two 
orders of magnitude larger than the other. 

Second, consider Solution 2, for which Eq. (|9b[) gives 
&s 1 s 2 — 0. (We have checked that even at the limits 
n\ = 80° and 100° of the allowed domain in Solution 
2, the proper spin-spin contribution is still one order 
of magnitude smaller than the self-interaction contribu- 
tion.) We conclude, that whenever one of the spins is 
aligned with the orbital angular momentum and the other 
spin is perpendicular to it, the proper spin-spin contribu- 
tion vanishes. This happens to be the case for Solution 
2 of the double pulsar J0737-3039. Thus <Jss-seif is the 
only spin-spin contribution in this case. We also remark 
that for Solution 2 the self-interaction spin contribution 



is only two orders of magnitude smaller than the spin- 
orbit contribution. This is because only the smaller spin 
contributes to (3. 

The importance of the self-interaction spin parameter 
is in gravitational wave astronomy, rather than in the tra- 
ditional electromagnetic observations. In ^(J the phasing 
function \1/ (/) was derived up to 2PN order. The coef- 
ficient a there denotes the proper spin-spin parameter. 
However as demonstrated here, the additional self-spin 
contributions given by crss-seif is of the same order of 
magnitude, thus it should be also included (together with 
gqm and <jdd, for high mass quadrupole and magnetic 
dipole moments, respectively). 

The self spin-spin contributions derived in this paper 
are more important whenever one of the spins is negligi- 
ble compared the other. To leading order in the (dom- 
inant) spin and when the mass of the high-spin com- 
ponent also dominates over its companion's, this is the 
Lense-Thirring approximation. Its first correction in the 
gravitational radiation is represented exactly by the SS- 
self contribution of the higher spin. 
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